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1 Introduction 



It has been noticed for a long time now, that the so called " Two-Hermitean-Matrix-Model" (in- 
troduced in particular for quantum gravity [1211101) and the so called " Complex-Matrix- Model" 
(used in particular for its applications to Laplacian growth models [2I1I2E], and string theory [I]) 
share lots of similarities: They have the same leading large N expansion properties, and, both are 
associated to some ensembles of biorthogonal polynomials which have formaly the same properties. 
Here, we add a new piece to make this correspondence more precise, we prove that both models 
have the same loop equations. 

Both models are not defined for the same weights, in fact, the set of weights for which one 
model is well defined has no intersection with the set of weights for which the other model is well 
defined. However, each model can be analyticaly continued to a larger set of weights, and in that 
sense, the two models coincide. 

When written in terms of eigenvalues, this identification of the 2-hermitean-matrix-model and 
complex-matrix-model has some interesting corolary: it gives a formula for computing integrals (of 
the Itzykzon-Zuber type) over the unitary group, as gaussian integrals over triangular matrices. 
Therefore, we obtain a very explicit formula for all correlators of the Shatashvili's type [21]. In 
[21] S. Shatashvili found a formula for all U(n) correlation functions, but his formula still contains 
integrals, is not explicitely symmetric in all variables, and is very difficult to use for practical 
purposes, such as 0. In the particular case of the 2-point correlation function, Morozov has 
found a much simpler formula [2B]. In j2B] A. Morozov computed it for U(n) with n < 3 and 
conjectured it for n > 3. Morozov's formula was later proven for all n in [Oj, and written in an 
even simpler form ^3]. Here, we find a natural generalization of Morozov's formula. The formula 
we find here, contains no integration, it gives the U(n) correlation functions as the sum of a finite 
number of terms, and is very efficient for effective computations. It also provides an alternative 
new proof of Itzykzon-Zuber's formula. 

The derivations proposed in this article are elementary, and it would be interesting to put 
them in the more general framework of group representation theory [201 E] • 

The main results presented in this paper are: 

• Theorem 13.31 and in particular Remark 13.31 which states the equivalence between the 
Hermitean-2-matrix model and the complex-matrix model: 

f dM 1 dM 2 F(M 1 ,M 2 )e-^ lTrMlM2 = [ dZ F(Z, Z ] ) e" 7Tr zz ' (1-1) 

JH n xH n JGL n (C) 

The definitions of each terms and the meaning of that equality are explained in section 13.31 

• Theorem 14. 1| which allows to compute U(n) integrals as triangular integrals. 

/ dUF(X,UYU^)e- TlXUYU] 

Ju(n) 
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E. E.(-ir(-l) T e-™ j T{n) F{X a + T,Y T + Tt) e~^t 
* A(X)A(Y) (1 " 2) 

for any polynomial invariant function F. 

• Theorem 15.21 which gives a formula for computing triangular matrix gaussian integrals. We 
parametrize polynomial invariant functions by pairs of permutations (of some size R), and a basis 
is written F^y. Theorem 15.21 gives the result of integration over triangular matrices: 

J T{n) dT e- Tr TT ' F n y g, y , X + T, Y + Tt) 

f JTp-Tr TTt 

JT(n) Ujt e 

= (^^(^y^n,^)^^,^^!,^^) ...^(^i/^l^l))^ (1-3) 

where M.( R '(x,y,X n , Y n ) is the matrix of size i?!, indexed by pairs of permutations: 



Theorem 15.31 shows that the matrices in eq 11-31 commute together, and can be simultaneousy 
diagonalized. 

• Theorem 16. 11 which gives a formula for computing correlation functions in terms of biorthog- 
onal polynomials: 



L „ rfM 1 rfM 2 F,^(f,y,M 1 ,M 2 )e- Tr ( yi ( Ml ) +y2 ( M2 ) +A//lM2 ) 

J Jti w, X Jti fi ,v 



(Mdet (MWfr&Q,!*)))^ 



L „ dM 1 dM 2 e- Tl ^( M ^ +v ^ M ^ M ^ 

J H n X ri n 

(1-5) 

where notations are explained in section 16.21 
Outline: 

• In part 2 we give definitions of groups and measures. 

• In part 3, we prove the equivalence between the Hermitean- 2- matrix model and the complex- 
matrix model, by showing that they have the same loop equations. 

• In part 4, we prove the identity between U(n) integrals and triangular integrals, and give 
some examples. In particular we rederive Itzykson-Zuber's formula and Morozov's formula. 

• In part 5, we compute the triangular integrals, by parametrizing polynomial invariant func- 
tions with pairs of permutations. In particular we explicit all four point functions. 

• In part 6, we integrate over eigenvalues using biorthogonal polynomials technics, and get 
expressions for correlation functions. 
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2 Definitions 

2.1 Ensembles 

Let 

• U(n) := group ofnxn unitary matrices, with the normalized Haar measure. 

• H n : = group ofnxn hermitean matrices, with the Lebesgue measure: 

dM := n dMu n dReM v dimM v 

i i<j 

• GL n (C) := group ofnxn complex matrices, with the Lebesgue measure: 

dZ := Y\ dKeZij dlmZij 

• T n :— group of n x n strictly upper triangular complex matrices, with the Lebesgue 

dT := JJ cZReTtf cflmT^ 

i<j 

• D n {R) := group of n x n real diagonal matrices, with the Lebesgue measure: 

dX =Y[dx u 

i 

• D n (C) := group of n x n complex diagonal matrices, with the Lebesgue measure: 

dX := Y[ dReXa dlmX u 

i 

• S(n) := group of permutations of n elements. 

2.2 Vandermonde determinant 

For any diagonal matrix X = diag(X 1; . . . , X n ) G D n (C), one writes: 

Apr) : n (A '' A - ! 

i<j 

and, for any permutation a G S(n), we define the diagonal matrix: 

X a := diag(X . (1 ), . . . , X a ^) 

Notice that: 

A(X a ) := (-If A(X) 



2.3 Invariant functions 

Definition 2.1 F(A,B) defined on GL n (C) x GL n {C) ^ C is an analytical invariant function 
if: 

• F is analytical in each variable, 

• W G GL* n (C), F(UAU-\ UBU- 1 ) = F(A, B). 

Examples: 

v / Rt \ 

F(A, B) = J] Tr J] (x t>rt - A)(y ttn - B) (2-9) 

t=l \n=l J 

F(AB) = e- M '^e- M2 f B > (2-10) 

Definition 2.2 Monomial invariant functions are functions of the form: 

P f Rt \ 

f(a,b) = n Tr n^ fct,rt Bh ' rt> > ( 2_n ) 

t=l \r t =l / 

where the k t)Tt 's and l t)Tt 's are integers such that k tjn + h,r t > 0- The total degree is 

p Rt 

degF:^J2J2 k ^ + kn (2-12) 

t=l r t =l 

Definition 2.3 Polynomial invariant functions are finite complex linear combinations of mono- 
mial invariant functions. 

Examples of polynomial invariant functions: 

F(A, B) = Tr A kl B h A k2 B h , F(A, B) = (l + Tr A kl B h ) (l + Tr A k ' 2 B h ) (2-13) 

p <i 
F(A, B) = JJdet(x t - A) kt l[det(y u - B) 1 - (2-14) 

t=l u=l 

F(A,B) =det(A® 1-1® B) (2-15) 

2.4 Decompositions 
2.4.1 Diagonalization 

It is a standard result in algebra (see [221 El 120] for instance), that any hermitean matrix M G H n 
can be written: 

M = UXU ] (2-16) 

where U G C/(n) and X G -D n (R). 



The measure is then: 

dM = J n A 2 (X)dUdX (2-17) 

where the Jacobian is 

n(n — 1) 
7T 2 

J n = , — (2-18) 

Uk=0 kl 

This decomposition is not unique. It is unique up to a permutation of eigenvalues, and up to 
multiplication of U by a diagonal matrix whose elements are on the unit circle. In other words, 
M = UXW provides a mapping between H n and U(n) x D n (R)/(U(l) n x S(n)). 

2.4.2 Jordanization 

A less standard result (see [2212312011111 for instance), is that any complex matrix Z e GL n (C) 
can be written: 

Z = U(X + T)U ] (2-19) 

where U E U(n), T G T n and X e D n (C). 
The measure is then: 

dZ = J n \A(X)\ 2 dUdTdX (2-20) 

where the Jacobian is 

n(n— 1) 

Jn = (2-21) 

-r-rn— 1 Ti \ / 

n fc =o k] 

This decomposition is not unique. It is unique up to a permutation of eigenvalues, and up to 
multiplication of U by a diagonal matrix whose elements are on the unit circle. In other words, 
Z = U(X + T)W provides a mapping between Gl n (C) and U(n) x T n x D n (C)/(U(l) n x E(n)). 

3 Gaussian matrix integrals 

In all what follows, we consider 3 complex numbers a±, a 2 and 7, and we define 

5 := a\Oi2 — 7 2 (3-1) 

and assume that 5^0. 

3.1 Gaussian Hermitean model 

Consider the measure on H n x H n : 

e _ Tr ^ M?+ ^ Mi+lMlM2) d ^ dm (3 _ 2) 
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Definition 3.1 The partition function is: 

Z ff (n, 7 ,Oi,a 2 ):= f dM, dM 2 e~ T "^ M ? + ^ M ^ M ^ (3-3) 

Notice that the integral Z H is absolutely convergent only if 

V0 E R Re(«ie + a 2 e" ± 2 7 ) > (3-4) 

which implies that Reai > 0, Rea: 2 > 0, (Re 7 ) 2 < Rea:iRea: 2 . 
An easy gaussian integral computation gives: 

,2 



n 



z " = r \n) ' (3 " 5) 

Definition 3.2 The expectation value of an invariant function F(A, B) is: 

J Il n X fin 

Remark 3.1 It is clear, from Wick's theorem, that if F is a monomial invariant function, then < F >h 
is a polynomial in ^ and j, and can be analiticaly continued to every complex 01,02,7, provided 
that 5/0. 

3.1.1 Gaussian Hermitean loop equations 

Consider a monomial matrix valued function, of the form: 

v Rt 
f(A,B) = f (A,B)l[Trf t (A,B), W = 0,...,p, f t (A, B) = J] A k ^B 1 ^ (3-7) 

t=l rt=l 

define: 

G (A, B) := J] Tr f u (A, B) , and if t > 1 , G t (A £) := J] Tr / U (A, £) (3-8) 

Theorem 3.1 One has the "loop equations": 

a, (G (Mi, M 2 ) Tr M 1 / (Mi, M 2 )) H + 7 (C (Mi, M 2 ) Tr M 2 / (Mi, M 2 ))) H 

Ro h),r~l I / /r-1 

lG G {M u M 2 ) Tr J]Mf°'"M 2 ^ ) Mf 

r=l j=0 \ V \u=l 

Tr (M 1 fco ' r ~ i ~ 1 M 2 io -'- ( Yl M^°' u M 2 l °-' 

\ \ u =r+l / / / jj 

p Rt kt,r~ 1 / / /r—1 

kt,r—j — l T./T It-, 



+ (G t (M u M 2 ) Tr m .U' .I// ) M(f (M 1 , M 2 )M^ r 3 X M 

t=l r =l j=0 \ \ \m=1 
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Rt 

fkt,u j, it h,u 



(3-9) 

and 



Yi .\/[" ■ .1 1 j 

\u=r+l / / I H 



(3 - 10) 



a 2 (G? (M 1 , M 2 ) Tr M 2 / (M 1; M 2 )) H + 7 (G (M 1; M 2 ) Tr M 1 / (M 1 , M 2 ))) H 

R «o,r-l / / /r-1 \ 

^ (g (M 1 ,M 2 ) Tr JjMf°' tt M 2 ^ M^°' r M 2 j 

r=l j=0 \ \ \u=l J 

Tr ( M 2 lo ' T ~ j ~ l ( JJ M^°' u M 2 l °- 

\ \u=r+l / / / H 

p R t h, r -l I / /r-1 \ 

+ \Gt{M x ,M 2 ) Tr m.\/?' .\/./ M 2 j f (M u M 2 )M 2 h ^- x 

t=l r=l j=0 \ \ \u=l J 

Rt 

Y[ M^ u M 2 l 



l . u 



Notice that the RHS is a linear combination of invariant polynomial functions of degree strictly 
lower than the LHS. 

Loop equations are a standard method for finding recursion relations among expectation values 
[HI] , they were first studied by [25] for the 2-matrix model, and solved more explicitely by [H3 H21 

Ej. 

proof: 

Write that the integral of a total derivative is zero: 

° = E/ dM ' dM "dk~. (^(M^M^e-^^^^ 2 ^^)) (3-11) 



i.e. 



— f..(Mi ( ,-Tv(^Mf+^-Ml+ 1 M 1 M 2 ) 



j dM,dM 2 (^— f^M u M 2 Yj e 
J dM 1 dM 2 f i>i {M 1 ,M 2 ) ( ai M Ui + 7 M 2ii ) e 



Tr (i^A^ + ^Mf + 7 M!M 2 ) 



(3-12) 

Similarly: 



I dM,dM 2 (^~^ (Ml ' M2) ) e~^^ M ^ M ^ M ^ 

i<j ^ u i ' 

J dM x dM 2 /<j(Mi, M 2 ) (ai(M^ + M Uj ) + j{M 2ji + M 2lj )) 



e _ Tr ( ^ Ml + Ml + T Mi M 2 ) 



(3 - 13) 



and 



I dM,dM 2 (gAf-.fiAMi,M 2 )) e -^(^ + ^M! + ,M lM2) 
£ / rfM i rfM 2 M 2) (ai(M ljt - M Uj ) + 7 (M 2ji - M 2y )) 



e _ Tr (^M 1 2 +^M|+ 7 MiM 2 ) 

(3 - 14) 

Taking + -i HEH , we get: 

i— — — - — ^ f ij (M 1 ,M 2 )] e -^^Ml + ^Mi +1 M lM2 ) 



dReM Uj dlmMuj 



+\J2J dM ± dM i 

+ i|/,M lrf M 2 ((^ + ^) /„<*,*)) e-Hr-MM*** » 

dMi<fM 2 (Tr/(M 1 ,M 2 )(a 1 M 1 + 7 M 2 )) e -Tr(^+ff m^m.m,) 



(3-15) 

i.e. one can proceed as if all the Miy were n 2 real indepedent variables, i.e., by abuse of notation 
we write: 



J2 f dM l dM 2 (y-^—f^M^M^ e 
dM 1 dM 2 ( Tr /(Mi, M 2 )( ai Mi + 7 M 2 )) e 



Tr (^M 2 +^Mf +7M1M2) 



Tr (^Af 1 2 + ^M|+ 7 M 1 M 2 ) 



(3 - 16) 

Now, one can use the following rules: 

• split rule: if /(Mi, M 2 ) = AMf 5 (where A and S are matrices), one has: 

E df{ sM^ h = E ■* (^f- 1 ") ^ M B ) < 3 - 17 ) 

i,j 1 7 1=0 

• merge rule: if /(Mi, M 2 ) = ATr (M^B) (where A and i? are matrices), one has: 

E = E Tr (^Mf-'-'BM!) (3-18) 



Then, if A and B depend on Mi, one has to use the chain rule. 

When one considers / given bv 13-71 one gets eq 13-91 

Eq J3-10I is obtained by doing the same for M 2 . □ 

We find again that (F) H is a polynomial in s& and ?. 

3.2 Gaussian Complex model 

Consider the measure on GL n {C): 

e -Tr^Z* + ^ +lZ *) dz (3 _ 19) 

Definition 3.3 The partition function is: 

Z (n,7,ai,a5i):= / dZe^ C^+^Vyzst) ( 3 -20) 

JGl n (C) 

Notice that the integral Zc is absolutely convergent only if 

V0 G R Re(«ie ie + a 2 e~ id + 2 7 ) > (3-21) 

One can see that with 9 = ir, this condition can never be compatible with 13-41 (with <fi = 0). 
Therefore, if Z# is an absolutely convergent integral then Zc is not, and vice- versa. 
An easy gaussian integration gives (where 5 = a\a 2 — 7 2 ): 

7T 



(3-22) 



which can be analiticaly continued to every ai,a 2 ,'~f, provided that 5^0. 
Definition 3.4 The expectation value of an invariant function F(A, B) is: 

f Gln(c) dZ F(Z, Zt) e" Tr (^^+^zt 2 +7 zzt) 
{F) ° " 2 J Me-K&^z^z*) (3 " 23) 



Remark 3.2 It is clear, from Wick's theorem, that if F is a monomial invariant function, then < F >c 

SSL 22 7 

«5 ' <5 ana s- 



is a polynomial in ^ and j, and can be analiticaly continued to every complex 01,02,7, provided 



that 5/0. 

3.2.1 Gaussian complex loop equations 

Consider a monomial matrix valued function, of the form: 

v Rt 
f(Z, Z+) = / (Z, Zt) J] Tr / t (Z, Zt) , / t (Z, Zt) = J] Z**-Zt^ (3-24) 

t=l r t =l 



define: 



G (Z, Zt) := J] Tr f u {Z, tf) , and if * > 1 , G t {Z, tf) := J] Tr f u (Z, tf) (3-25) 
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Theorem 3.2 One has the same loop equations than theorem \S.l\ with replacing the subscript H 
byC. 

a x (G {Z, tf) Tr Zf {Z, tf)) c + 7 {G (Z, tf) Tr Z*f (Z, tf))) c 

R fco,r-l / / /r-1 \ 

r=l j=0 \ \ \u=l / 



Tr \ Z h °--i- x Z^' r ( f[ Z k ^Z^ l0 Aj\ 



c 



V Rt h,r-l I / /r-1 

+ \ G t {Z,tf) Tr m^'^Zt^ ) ZifoiZ^Z^Z^-^Z^ 

t=l r=l j=0 \ \ \«=1 



(3 - 26) 

and 



v u=r+l / / 



r* 



a 2 <G (^, Z f ) Tr Zt/ (Z, ^)> c + 7 <G (Z, ^) Tr Z/ (Z, Z+))> c 

R lo,r-l I I /r-1 

( G o(Z,Z^) Tr Y\Z k ^Z ]khu ) Z ko ''ZV 

r=l j=0 \ \ \u=l 

Ro 

7 j;l0,r-j-l TT 7kou7j; lo,u 



Tr Z^' 3 ' 1 11 Z k ^Z^ 



(3 - 27) 



\,M=r+l J J I (j 

V R t h,r-l I / /r-1 \ 

t=l r=l j=0 \ \ \U=1 

Rt 

\M=r+l / / I Q 



h,r-j-l 



Notice that the RHS is a linear combination of invariant polynomial functions of degree strictly 

lower than the LHS. 

proof: 

The proof is very similar to that of theorem 13.11 Write that the integral of a total derivative 
is zero: 

(3 - 28) 



and 
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(3 - 29) 



Taking the sum of both lines, one can proceed as if all the Zij and Z\j were real indepedent 
variables, and from there, follow the proof of theorem 13.11 □ 

Remark 3.3 We see that the loop equations of both models are identical. It is clear from the above 
derivation that this is general, even for non gaussian measures. When the measure is gaussian, the 
loop equations determine completely every expectation value, while for non-gaussian measures, the loop 
equations give recursion relations for expectation values, but don't give the initial conditions. 
Let us consider in particular the "semi-classical case" jUEj, i.e. with a measure of the type 

dfi(Mi, M 2 ) = e~ Tr [v 1 (M 1 )+y 2 (M 2 )+A/ 1 Af 2 ] (3 . 30) 

where V[ and V^' are rational functions. In that case, the initial conditions which allow to determine all 
polynomial expectation values recursively, are in one-to-one correspondance with homology classes of 
integration paths for pairs of eigenvalues [7j, therefore, there exists a choice of integration path T such 
that one can write: 

/ dMi dM 2 e" Tr [Vi(M 1 )+v 2 (M 2 )+Af 1 M 2 ] _ f dz Q ~ Tr [l/ 1 (Z)+V 2 (Zt)+ZZt] £g_ 31 j 

J{H n xH n )(T) ' JGL n {C) 

and one can consider that this equality defines the RHS. Somehow, the complex matrix model is nothing 
but the analytical continuation of the 2-matrix model defined on some classes of contours. 



3.3 Relation between the two models 

Theorem 3.3 For any polynomial invariant function F(A,B), one has: 

(F) H = (F) c (3-32) 

Notice that (F) H and (F) c have been defined for different range of values of a±, a 2 , and 7, 
but, as we have explained above, both are polynomials of 2& and % (and can be analyticaly 
continued to any a>i, a 2 , and 7). Theorem 13.31 is thus an equality between polynomials, 
proof: 

It is sufficient to prove it for monomial invariant functions. The proof is clearly obtained from 
the loop equations, by recursion on degF. It is obviously true for degF = 0, i.e. F — 1. And the 
loop equations of both models are identical. □ 

Definition 3.5 For any two given complex diagonal matrices X and Y , and any polynomial 
invariant function F, define: 



W F (X,Y) := A 2 (X)A 2 (F) / dU F(X,UYU^e-^ XUYU ^ (3-33) 

JU(n) 

L< , dT F(X + T,Y + Tt) e -^ TrTTt 

u F (X, Y) := A(X) A(y) ^ , Tp -^TTt ( 3 - 34 ) 

JT(n) al ^ 
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which is a polynomial in all its variables Xi, Yj, and a polynomial in and: 
iy F (X,y):=^VVA(X (T )A(F T )e-^ TrM f dT F(X a + T,Y T + Tt) e ^ TrTTt (3-35) 

U - a T JT(n) 

Theorem 3.4 For any polynomial invariant function F(A,B), one has: 

, 2 = /" r / dXdYe-^ x2 e-^ Y2 W F (X,Y) 

n\ 2 Z H {n,^,a 1 ,a 2 ) J Dn (K)xD n (R) 

= -^rj^ r / dXe-^ x2 e-^ 2 W F (X,X) (3-36) 

n\Z c {n, 7, ati,a 2 ) J Dn {c) 

proof: 

Start from theorem l3.3[ diagonalize M\ and M 2 on the hermitean side, and jordanize Z on the 
complex side. 

(F) H = — /" T / dXdY e-^ TrX2 e-^ TrY2 W F (X } Y) 

t (R)x»„(rt) 

(F) c = yry ( ° n / dX e-^ x2 e-^ T ^\-^ xx uj f (X,X) 

*(C) 



p 


n\ 


2 Zij(n,7, ai,a 2 ) 




Jn 


n\ 


Zc(n,7,ai,a 2 ) 






n\ 


Z c (n,"/,a 1 ,a 2 ) 




Jn 


n\ 


Zc(n,j,ai,a 2 ) 



T / dX e-^ x2 e-^ Trx2 e^ TTX ^uj F (X„X T 

JD„(C) 

[ dX e-^ x2 e-^ x2 W F (X a} X T ) 

Jd„(C) 



(3 - 37) 

The equality in the first line is obtained by diagonalizing Mi and M 2 (with Jacobian given in 
eq 12- 17(1 . the equality in the second line is obtained by Jordanizing Z (with Jacobian given in 
cq l2-19j) . the equality between the second and third line holds for any pair of permutations a and 
r (it can be proven with the Lemma lA.ll given in appendix), and the equality of the last line 
comes from the definition of W F . □ 

4 Unitary group integrals 

Here is one of the most important theorems of this paper: 

4.1 Unitary integrals and triangular integrals 

Theorem 4.1 For any invariant function F(A,B) one has: 

[ dUF(X,UYU^)e-^ XUYW 

JU(n) 

c n E T (-ir(-l) T e-^ M J T(n) F(X a + T,Y T + Tt) e"^^ dT 



n\ A(X)A(Y) 
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where 



(4-2) 



i.e. 

W F (X, Y ) = n\ Cn W F (X, Y) (4-3) 

proof: 

Using the Lemma lA.ll given in appendix, and using theorem 13.41 we have: 
/ dXdY e-^ x2 e-^ Y2 W F (X,Y) 

/ dX e-^ TrX2 e-^ T ^ 2 W F (X,X) 

JD n (C) 

±=Y [ dX e-^^ x2 e-^ T ^e-^ x ^MX„,X T ) 



Zh Jn 



Jl n\ Z c 



-^% J "^ / dXdY e-^ x2 e-^ Y2 e-^ x ^MX„Y T ) 

^ ffi I dXdY e-^ x2 e-^ Y2 W F (X,Y) 

2 2*\ n\Z C JD n (R)xD n {R) 



(4-4) 



J2 (i*J J D n (R) x D„ (R) 

n\c n [ dXdY e-^ TrX2 e-^ Try V F (X,F) 

J D„(R)xD„(R) 



Notice that if /(A) and g(-B) are invariant functions i.e. f(UAU x ) = f(A) for all A and {7 (resp. 
giUBU- 1 ) = g(B) for all 5 and U), one has: 

W) (X Y) = f(X)g{Y)W F (X, Y) , W fix)g{Y)nX)Y) (X, Y) = f{X)g(Y)W F (X, Y) 

(4-5) 

Thus, for any / and g: 

0= f dXdYe~^ i:tx2 e-^ TTY2 f(X)g(Y)(n\c n W F (X,Y)-W F (X,Y)) (4-6) 

J D„(R)xD„(R) 

Since W^p(X, K) and n\c n W F (X , Y) are symmetric functions and entire functions of all their vari- 
ables, they must be identicaly equal to one another. □ 

4.2 Examples 

Let us illustrate theorem 14.11 on some simple examples and recover some classical results. 
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4.2.1 Harish-Chandra-Itzykson— Zuber's formula 

We can use theorem 14.11 to find a new proof of the famous Harish-Chandra-Itzykson-Zuber's 
formula [X31 12D| - 

Indeed, consider F(A,B) = 1, theorem 14. II gives: 

n V ( —~\ ) a ( —~\ ) T TT p~~i x °i Y n r 

->yTr XUYW C ™ ^cr,r\ X J V ) 1 Li c / , rp --y Tr TT+ 

un _ " n! A(X)A(F) ' " r 

n ( n — 1 ) 

tA^ - det£ 



''" A(X)A(F) 
(4-7) 

which is the famous Harish Chandra-Itzykzon-Zuber integral. Here E is the matrix 

Eij := e" 7 ^ (4-8) 

4.2.2 Morozov's formula 

Consider Tr A k B l for any integers k and I. It is in fact simpler to introduce a generating function: 

oo oo 

y k=0 1=0 y 

which is to be understood as a formal power series in its large x and large y expansion. F(A, B) 
is merely a convenient way of considering all polynomial invariant functions of type Tr A k B l at 
once. 

We have: 

x-(X + T) = ]^ {x~^X T ) V^X (4 " 10) 

and thus: 



f T{n) dTe-^TTt J T(n) x _ (X + T) y _ {Y + T t) 

n n p+1 ^ 9+1 ^ 

p=0 9=0 ii<i 2 ...<i p+ i ji<h-<j q +i k=l lk 1=1 y n 

JT(n) ± n,i2 ± 12,13 ■ ■ ■ ± ip,i P +l J -j q+1 ,j q • • • ± j 2 ,jl _ x 

/T(n)^-^TTt " ^ 

That last integral is non vanishing only if p = q, and according to Wick's theorem, it is the sum 
of all possible pairings. Because of the ordering of the z^'s and jVs, the only non vanishing pairing 
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is obtained for %k = jk for all k. Therefore: 



' f Tr 1 „„ dTe~^ TTl 



J T{n) dTe-^TTi J T(n) x _ (X + T) y -(Y + T^ 

oo 1 p+1 

ee in 

p=0 ii<i2— <i 



n 



-7 



i=l 



and then theorem 14.11 gives: 

n(n — 1) « / 1 1 \ 

(T) 2 / dC/Tr ( — —U— —rf ) e -^ XUYUi 
w J u(n) \x-X y-Y J 

n\ A(X)A(Y) 
- det E + det ( E + ± 



7Cn A(X)A(E) 

? (- 1 + det ( i + ^ ^^vi) - sot (4 - i3> 



i.e. 



( ) K " } -1 f-l + detfl + i^— E^—ET 1 ) ) (4-14) 



which is identical (for 7 = —1) to what was found in 0EI]> i.e. the compact version of Morozov's 
formula 



5 Computation of triangular integrals 

The goal of this section is to compute the triangular integral on the RHS of theorem 14.11 Here, 
we consider 7 = 1. 

5.1 Parametrization of polynomial invariant functions 

Definition 5.1 Let R be a positive integer. Let x = (x±, . . . , xr) and y = (y±, . . . , yn) be 2R 
complex numbers. Let tc and it' be two permutations o/E#. 

The permutation tttt'^ 1 is made of p cycles C±, . . . , C p of lenght . . . , R p which we note: 

C k = (ik,i^jk,i ~* k,2^jk,2 ^ ik,3~ • • ~* -^jk,R k ik,i) (5-1) 
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We define, for (A, B) G GL n {C) 2 in any dimension n: 



P / Rk 1 \ 

F^'(x,y,A,B) : = J] ^ + Tr JJ - _^ _ - (5-2) 

As explained above, this definition is to be understood as a formal power series in the large Xi and 
yj expansions, it is merely a way of considering all polynomial invariant functions at once. 

Examples: with R = 2, we have: 

F(i){2),(i){2)(x 1 ,x 2 ,yi,V2, A,B) = ( 1 + IV - 1 - - 1 - ) (l+Tr 



x 1 - Ay 1 - B J \ x 2 -Ay 2 -B 
^(12), (12) (xi, x 2 , yi, y 2 , A, B) = (l + Tr ^~^) ^ x '_ Ayi \ B ) 

F (1)(2) , (12) ( Xl ,x 2 , y u y 2 , A, B) = Tr Xi \ Ayi \ B X2 \ Ay^- B 
F (12) , (1)(2) ( Xl ,x 2 , y u y 2 , A, B) = Tr ^J^^J^ 

(5-3) 

Definition 5.2 Let R be a positive integer, x — (x±, . . . , xr) and y = (y±, . . . , yn) be 2R complex 
numbers. Let n and ir' be two permutations ofT> R . Let n be an integer, and X = diag(Xx, . . . , X n ) 
and Y = diag(Y"i, . . . , Y n ) be two complex diagonal matrices of size n, We define: 

Wfy (x, y, X, Y) := 1 if n = or R = (5-4) 

W^,(x, y, X, Y) := F n y(x, y, X u Y ± ) if n = 1 (5-5) 



and otherwise 



y, X, Y) : = ^ , ' (5-6) 



f TM dT e" Tr TT] F ny (x, y, X + T, Y + 7*) 



It(u) dT 



C" 



#ere ; j;-(x+t) ^ S defined by: 



U-i) 



f 1 \ ^ ^ ^ 1 J J 1 

^-(x+T)J i7 . :_ + 2^ ^x: Tmi ^x~ Tji ' 12 • • • t^^t - v~ 
n ^ ' ' *j p=i t<ii<...<i P <j 1 p j 

(5-7) 

5.2 Computation of triangular integrals of invariant functions 

We are now going to find some recursion relation in n for the W's. 
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Theorem 5.1 



W*J (x, y, x,y)=y, m S ( f > y> **> Y ^ w pT } & y> x> Y ) 



(5-? 



where X := diag(X 1; . . . , X n _i), F := diag(F 1? . . . , F n _i) ; and: 



V, X n, Yn) = J] ^W.PW + 



1=1 



(xj - X n )(y 7r( j) - F n 



(5-9) 



proof: 

If T is a strictly upper triangular matrix of size n, we define T the triangular matrix of size 
n — 1, such that Tj j = 7^ for all z, j < n, and w the vector made of the last column of T, = T^ n : 



/•-. 



Ml 



Wn-1 
/ 



(5-10) 



We define 

Notice that: 



1 



3i,n&j,n 

x -(X + T)J tj \ x -(X + f)J id ' ^-Xn^\x-(X + f)J itk " K ' x-X n 



n-1 

£ 



(5-11) 



and 



1 + Tr 



1 + Tr 



1 



1 



x - (X + T) y - (F + Tt) 
1 1 

x - (X + f) y - (F + ft) 

/ n— 1 Ti—l 



(x-x n )(j/-Y n ; 



(5 - 12) 



rtr \y-(F+Tt)x-(x+T) 



Now, we integrate u out, using Wick's theorem, i.e. take the sum over all possible pairings of a u 
and a w. The pairing (uk,ui) gives a factor 5^. 

Let us represent Wasa bivalent graph G, whose edges are pairs (xi, y n (i)), and whose vertices 
are pairs {y^^^Xi). 

Relation eq 15-1 II means that, for each edge (xi,y n ^) of G, we can either: 

- let the edge untouched (first term in eq i5-ll|) . with weight 1, 

- remove the edge (second term in eq l5-llj) . with weight ( x ._ Xn )( y n -Y n ) > 
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- remove the vertex (y^ni), xA (third term in ea 15-1 ID . with weight -, tft, which 

means that either the neighboring edges cannot stay untouched. 

Then, we integrate u out, i.e. we take the sum over all possible pairings, i.e. we draw new 
edges between vertices (those not removed), so that the final graph is bivalent. For each pairing, 
we get a new graph G'. The sum over possible pairings, is thus the sum over bivalent graphs G', 
whose vertices form a subset of the vertices of G, i.e. 

= M G,G' W$~ x) (5-13) 

G' 

where the coefficient M.g,g' is computed as follows: 

- M.Q G' receives a factor 1 + -, ^ for each edge (xj, y^u)) of G which is unchanged, 

i.e. which is an edge of G'. (1 if it was not removed, and -, ^ if it was removed and 

drawn again). 

- the weight of each edge (xi,y n (i)) of G, which is not an edge of G', is ^ x ._ x ^ y _ Y s • 

- the weight of removing a vertex is the same as the weight of creating a lenght 1 cycle at that 
vertex. In other words, if G' has less vertices than G, consider G" obtained from G' by adding 
lenght 1 cycles at each missing vertex, one has M.g,G' — M-g,g>'- The sum over G' can thus be 
written as a sum over G" , where G" has as many vertices as G, and all cycles of lenght 1 come 
together with a 1 added. 

- relation eq J5- 121 ensures that the previous rules apply also when G has lenght 1 cycles. 
To sumarize, we have: 

w g ) = ^Mg,g» W$T x) (5-14) 



G" 



where 



Mg,g " n ( i+ u x n )( yMi) -Y n )) n u _ Xn )(y Mi) -Y n ) (5_i5) 

when G and G" are written in terms of pairs of permutations, it reduces to eq 15-91 □ 
Remark 5.1 Notice that: 

M^(x,y,X n ,Y n )=M^ R \x,y,X n ,Y n ) t (5-16) 

Mi R l(x, y, X n , Y n ) = MfX y ^ (y, x, Y n , X n ) (5-17) 

M^jpix, y, X n , Y n ) = M { *1 (£ p -i ,y,X n , Y n ) (5-18) 

Theorem 5.2 

W^,(x,y,X,Y) = (M {R) (x,y,X n ,Y n )M {R \x,y,X n . u Y n ^) . . .M {R \x,y,X 1 ,Y 1 ))^ l (5-19) 



18 



proof: 

For n — 1, we have 

W^(x, y, X u Yt) = F ny (x, y, X u Y 1 ) = M^x, y, X h Y x ) (5-20) 
The proof follows from recursion on n. □ 

Theorem 5.3 The matrices J^A^ R \x,y,^,rj) commute among themselves: 

M( R \x, y, £, V )M( R \x, y, rf) = M^ R \x, y, rf)M^(x, y, £, rj) (5-21) 

proof: 

Let n = 2, X = diag(-Xi, X2) and Y = diag(Y"i,Y2) be two diagonal matrices, and X = 
diag(X 2 , Xi) and Y = diag(Y2> Yi). Let T be a 2 x 2 upper triangular matrix with non vanishing 
element Let U be the 2x2 matrix: 

C=fv Tl \, Y '~ Yl ) (5-22) 



Xi — X2 T12 

it satisfies: 

U(X + T) = (X + T)U , U(Y + T f ) = (Y + T l )U (5-23) 
If U is invertible (which is true for almost every T), one has: 

F VtP (x, y,X + T,Y + T^) = F^ p (x, y,X + T,Y + T l ) (5-24) 



for every T (except a zero measure subset). Since the Jacobian 



c)T 
9T 



1, one has: 



J T{2) dTe-^ TT ' F^ p (x,y,X + T,Y + T^) _ J T(2) df e~^ ff1 F^ p (x,y,X + T,Y + T*) 

Using theorem 15.21 for n = 2, we have: 

M( R \x, y, X x , Y 1 )M^ R \x, y, X 2 , Y 2 ) = M {R \x, y, X 2 , Y 2 )M^(x, y, X x , Y x ) (5-26) 

□ 

Corollary 5.1 therefore, there exists an orthogonal matrix U(x,y) , independent of £ andr], such 
that: 

A(S,y } tv) :=U{x,y)M {R \x,y^,r,)U\x,y) (5-27) 

is a diagonal matrix 

A(x, y, f , 77) = diag (A 7r (f , y, £, 77)) (5-28) 
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Notice that A is a rational function of £ and r\. 
Thus: 

L, , dTe~ TrTT] F n n ,(x, y,X + T,Y + T^) ^ JL 

— _ TrTTt = £ U n , p (x, y) U^^x, y) J] A p (x , y, X t , Y,) (5-29) 

JT{n) a± e p i= i 

and: 

/,(„) ^ e- 1 *™» F„„ (f , g, X, UYW) Aet(e-x.r.AMy.X,Y l )) 
, i4 , e -«™. " det(e-W) 

Jt/(n) p V / 

(5-30) 

Remark 5.2 if one defines the "Matricial determinant" as follows: 

Definition 5.3 Let M G GL n {Gl m (G)) , i.e. for each i = 1, . . . , n, j = 1, . . . , n, Mij is a square matrices 
of size m. We define: 

1 n 
Mdet(M):=- £ £ (-1)^-1)^^(0,^) ( 5 ~ 31 ) 

o-eS(n) reS(n) i=l 

which is a m x m square matrix. 
then we have: 

n(n-i) ( Mdet (e' x ^ M ^ (x, y,X t ,Yj))) n 

au r n y{x,y,A,u i U')e = c n (ir; 

'U(n) 



f + -r Y rr V rrt n(n-i) ( Mdet(e _XlY ? y, Xj, / 

/ dUF n y(x,y,X,UYU^e-^ XUYUt = c n (vr)^ ^ ^ AnnXm ^ 



if it! = 0, one immediately recovers the Itzykson-Zuber's formula, and if R = 1, one immediately recovers 
Morozov's formula. 

5.3 Examples 

• Example R—l: 

M$(x,v,t,v) = l + X 1 _z y 1 _ r} (5-33) 



and thus: 



f ciT e~ Tr TTt (l + Tr * * ^ 

JT{n) al e ^± -T ±1 ^-(x+T) 2!-(Y+Tt) ^ 



+ (5-34) 



• Example R= 2: 
We have: 



%)(2),(i)(2)(zi, a?2, 3/1,2/2, A -B) = ( 1 + Tr - 1 A ^ 1 — J (l+Tr 



x 1 -Ay 1 -BJ \ x 2 -Ay 2 -B 
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F(i)(2),(i2)(xi,X2,y 1 ,y 2 ,A,B) = Tr 



xi - A y 1 - B x 2 - A y 2 - B 

F 0*UW)(W* VI. V* A B) = Tr ^Ay^-Bx^-Ay^-B 
(5 - 35) 

and 

- A ^a)(2),(l)(2) (^1^2, 2/1, 2/2, ^7) = + (l + ^^) 

^), (12) (^i^2,yi,y2,e^)= + + (5 _ 36) 

• M (i)(2),(i2) (^i, ^2, 2/i, 2/2, ^, ^7) = x ^ yi \ X2 ^ y2 \ 

, - A ^(l2),(l)(2)( a; i' a: 2,2/l,2/2,^,^) = X *_S x ^ yl_ r] 

i.e., the matrix M^ 2 \xi, x 2) y± : y 2 , £, rj) is: 



\ xi-£vi-vj \ vi-KVi-r) ) 



1111 



xi-£ Vi— i) x 2-t, Vi-r\ 

1 1 



(5-37) 



xi J/2-»7 a;2-C 2/1—1? 

i.e. 

1 / I .1 \ ( 1 . . / 1 

1 



2/i — ^ 2/2 - ^ 
1 /l + S" 1 



(xi -0(^2- 0(2/i -^)(y 2 - v) V 1 1-5 

(5 - 38) 

where 

S = ^i-x 2 )( yi -y 2 ) (5-39) 
Define the following orthogonal matrix (JA^ 2 \xi, x 2) y± : y 2 ) U^ 2 \xi, x 2) y±, y 2 ) 1 = 1): 

ZY (2) (zi,z 2 ,2/i,2/ 2 ) := -j^=== (^ s \ x X ~ 1 S ^) , whereA = v / TT^ (5-40) 

one has: 

M (2 \x u x 2 , y b j/2, »7) = U {2 \x u x 2 , y 1: y 2 ) A {2 \x u x 2 , y 1} y 2 , f , 77) U (2 \x 1: x 2 , y u y 2 f (5-41) 
where A (2) (xi, x 2 , 2/1, 2/2, v) = ^s(^+(x 1 ,x 2: y 1: y 2 ^,rj) : A_(x 1: x 2: y 1 ,y 2: ^rj) ) with 

A±(x 1 ,x 2 ,y 1 ,y 2 ,^,r]) = 1 + ^ ( — ^— r + — ^— r J ( — - — + 



2\xi-t x2-tj\yi- r l V2-r] / 
^(xi -0(^2 -0(2/i (2/2 -?7) ^ ^ 



Eventualy, one gets: 

JT(n) U1 ^ 11 X!-(X +T) ^-(y+Tt) s 2 -(X+T) y 2 -(Y+T1) 

fr(n) 
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f r/T e" Tr TTt 



1 f n n \ 

2A ( II A ~( x u x ^y^y2,X i ,Y i ) - Yl^+(xi,x 2 ,y 1 ,y 2 ,X i , Y$ J (5-43) 

\i=l i=l J 



/r(n) dTe i 1 + Tr x x -{X+T) yi-(Y+Tt)) + Tr X2 -(l 



(X+T) y 2 -(Y+Tt) 



JT(n) al C 

1 ( n H \ 

= 2x[(\ + S)Y[A + ( Xl , x 2 , Vl ,y 2 , Xi, Yi) + (A - S) JJ A- (an, x 2 , Vl ,y 2 , X h Y t ) (5-44) 

\ i=\ i=l / 

6 Mixed correlation functions and biorthogonal polyno- 
mials 

Let us consider two polynomial potentials V\{x) and V 2 (y) . Our goal is to compute the following 
matrix expectation values: 

j HnXHn dM 1 dM 2 F ny (x, y, M u M 2 ) e" Tr (Mmq+vmm^+a/^) 

JZ Z rfM 1 rfM 2 e- Tr W( M i)+^(M 2 )+M 1 M 2 ) ( 6 - X ) 

6.1 Biort honor mal polynomials 

We recall here a few elementary notions about biorthogonal polynomials. More detailed descrip- 
tions can be found in particular in [221 El El El El E] • 

We introduce two families of polynomials p n (x) = -k=x n + 0(x n ~ 1 ), q n (y) = ~k=y n + Oiy 12 " 1 ), 

V fin V^n 

with the same leading coefficient -4=, and orthonormal with respect to the pairing: 




(Pn,q m ) = / / dxdyp n (x)q m (y)e- {Vlix)+V2(y)+xy) = 6 nm (6-2) 



The integration path is a priori R x R, but this condition can be relaxed (see [HE]). When they 
exist, the biorthonormal polynomials are uniquely determined. 

Since the biorthonormal polynomials form a basis, one can decompose xp n (x) onto the basis 
of p m (x) with m < n + 1: 

n+l 

xp n {x) = ^Q nm p m {x) (6-3) 

m=0 

and similarly: 

n+l 
m=0 

Q and P are infinite matrices. In the case where V 2 (resp. Vi) is a polynomial, then Q (resp. P) 
is a finite band matrix. 



22 



We also introduce the following ooxn rectangular matrix: 

/I \ 



n n _i := 



1 



V 



(6-5) 



which is the projector onto the n first polynomials. 

6.2 Mixed correlation functions 
Theorem 6.1 

f H xH dM 1 dM 2 ^ x (f,y,M 1 ,M 2 )e- Tr ^( Ml ) +V2 ( M2 ) +MlM2 ) 



L „ c /Mic/M 2 e- Tt W( M i)+V2(M 2 )+M 1 M 2 ) 

J ii n X H n 

U*, P & V) ^n', P {x, y) det (n^! : A p (x, y, Q, P*) : U n ^) 



(6-6) 



where for any function of two variables f (£,,?]), we define : f{Q ) P t ) : by putting the Q's on the 
right of the P's. This is always possible in this case because A p (x,y,£,r)) is a rational function of 
£ and rj. 

proof: 

It works as usual (see |2*T) 122]). by writings Vandermonde determinants as: 

n-l 

A(X) = detiXf 1 ) = det(^h~ lPj ^(X t )) = l[Vh~ ]T(-ir n^MTO ( 6 - 7 ) 



i=o 



A(Y) = detiYr 1 ) = det(VvTg i _ 1 (y i )) = J] V* ^(-1)^^)^) 

j=0 t i 

Then, we use eq J5-3Ul i.e. 

dM 1 dM 2 F ny (x, y, M u M 2 ) e~ Tr (WHWHm^) 



(6-8) 



~ n(n-l) . 

n-1 



1 " x 

— JI^ E ^A^y)^'A^y) E 

n— 1 n 

— n ^ e - ^ f > $ e n : ' & q> pt ) 



=0 pes(ij) 



<T,T,i/6S(n) 



i=l 
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n— 1 n 

H hi Yl Z4,pOE,y)Z4',p(z,y) J] ( _1 ) CT II : A p&V'Qi pt ) ; i,<7(i) 

j=0 peS(_R) creS(n) t=l 

n-1 

II ^ 5Z ' $ $ det (X-l : A p( f ' £ Q> Pt ) '■ n "-l) 

(6-9) 

□ 

or, using the matricial determinant defined in def J5.3l 

J HnXHn dM x dM 2 F wy (x, y, Mi, M 2 ) e" Tr (WHWHM^) 

f„ ~ dM l dM 2 e-^W (M 1 )+V 2 (M 2 )+M 1 M 2 ) 

= (Mde^(Ul. 1 :M^{S,y,Q,P t ):U n -i)) VtVl (6-10) 

Example: with i? = 1, we find: 

f„ „ dMi dM 2 (1 + Tr ^rr^rr) e-^^^H^^J+JifiJWa) 
f„ ~ dM 1 dM 2 e- Tt ( v ^ M ^+ v ^ M ^+ M ^ 

which is identical to what was found in [Hj. 

7 Conclusions 

In this article, we have shown that the hermitean 2-matrix model and the complex matrix model 
have the same loop equations. In the gaussian case, that implies they are identical. In case the 
weight is non-gaussian, the loop equations, which are recursion equations, determine all correlation 
functions when some initial conditions (moduli) are fixed. The generalization of the hermitean 
2-matrix model to homology classes of contours (as in [S]), allows to have any arbitrary initial 
condidtions, so, there exists a choice of homology class of contours for each set of initial conditions, 
i.e. for which the complex matrix model is identical to the 2-hermitean matrix model. Conversely, 
the initial conditions for the complex matrix model are not fully understood yet, they depend on 
how the complex matrix model is defined. If the complex matrix model is only a formal integral 
defined by its large n properties as in j23|2Hj, initial conditions are associated to filling fractions, 
and can thus be chosen arbitrarily. If the complex matrix model is defined as the result of a 
convergent integral for all n, it is not known yet how to find which homology class of contours it 
corresponds to. 

The consequence of that identification, through diagonalization of hermitean matrices and Jor- 
danization of complex matrices, yields an identity between unitary group integrals and triangular 
matrices integrals, which seems to be a special case of the identification of GL n (C)/T(n) and 
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the quotient of SU(n) by its Cartan subalgebra. The nature of that identification needs to be 
further understood, in particular in terms of characters of both groups, and in terms of group 
representation theory, in terms of Weyl's character formula, or Harish- Chandra formulae. 

The gaussian triangular matrix integrals are easily computed, and we thus get very explicit 
expressions for all expectation values of the type which were studied by Shatashvili [21]. In 
particular, we have provided a new proof of the Itzykson-Zuber-Harish-Chandra integral, as well 
as Morozov's integral. The key piece in this computation is that the matrices M. commute 
together. This fact seems to be related to some Yang-Baxter relations, and it would be interesting 
to understand how. 

It would be interesting also to understand these formulae in the framework of Duistermaat- 
Heckman semiclassical theories [TT] . 

Then, we have been able to perform the integral over eigenvalues, in a way very similar to 
what was done in [B], i.e. in terms of n x n determinants. It would then be interesting to rewrite 
these n x n determinants in terms of determinants of size independent of n, using kernels, as it is 
known for non- mixed expectations values (see 0IH1EII])- 
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Appendix A Gaussian integrals 



Let 6 = aia 2 — 7 2 . 
Real integrals: 



/ dxdye-^ 2 ^y' 2+ ^ = ^ (1-1) 



/rxr dx dy x k y l e ( 2 X + 2 v +i x v) 
J RxR dxdye~^ x2+ ^ y2+ ^ 



if k + lis odd 




k-l 



da 2 J V ®1 J \VS 



L 



D„(R)xD„(R) 



dXdYe- Tl ^ x2+ ^ Y2+ ^ = (^=J (1-3) 
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Complex integrals: 



dxe~^ x2+ ^ x2+ ^ = —— (1-4) 



/*, dx x k x l e ( 2 x + 2 X +i xx ) 

' = if k + 1 is odd 



ot-\ 9 , an n 

— L^J — ^* | ,., ., , 

fc-Z 



j'dxe ^2 x2 + 2 1 +7^) 



, «— i 7 

8\ f „ d \~ { d\ l 



2— -— -= if k > I 
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l-k 



5\ f „ d \ ~ f d 



2— -— -== if k < I 



7T 



7T / \ da 2 J V ^7 / V \/— $ 
(1-5) 

a^e"^^^^ = f -£=Y (1-6) 



D„(C) \V— 0/ 

Lemma A.l Lei y), fre a polynomial in all its variables X\, . . . , X n and Yi, . . . , y n; one 
has: 

(1-7) 

proof: 

Eqs 11-21 and 11-51 show that it is true for n — 1. By decomposing w into monomials, the integral 
decouples into a product of n = 1 type integrals. □ 

Appendix B Some Commutations 

Theorem B.l Tne matrix A4^(x,y,^,r]) commutes with the matrix A(x,y) defined by: 

Air^yX, if) . Xiy^a) 

A^y(x, y) '■= 1 if 7T7r /_1 = transposition (2-1) 
A^y(x,y) := otherwise 

Theorem B.2 The matrices A a,l3 (x,y) defined by: 

A"'P,(x,y) := Sp y a ) TT [S^y^ H J 

it V x a -xiy p - y 



1-5 



{x ct -x^i w )(y p -y <ct) ) . s ^ , 
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][ ( <^7r(i),7r'(i) H ) 



(2-2) 



commute together for all a, f3. They also commute with AA(x, y, £, 77) and with A(x, y). 
One has: 

M (R \x , y, £, 77) = 1 + 77 1 7 ^(S, y) (2-3) 
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